On the stability of spherically symmetric spacetimes in metric f(R) gravity 
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We consider stability properties of spherically symmetric spacetimes of stars in metric f(R) grav- 
ity. We stress that these not only depend on the particular model, but also on the specific physical 
configuration. Typically configurations giving the desired 7ppn ~ 1 are strongly constrained, while 
those corresponding to 7ppn ~ 1/2 are less affected. Furthermore, even when the former are found 
strictly stable in time, the domain of acceptable static spherical solutions typically shrinks to a 
point in the phase space. Unless a physical reason to prefer such a particular configuration can be 
found, this poses a naturalness problem for the currently known metric f(R) models for accelerating 
expansion of the Universe. 
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I. INTRODUCTION 



The observation that the expansion of the Universe ap- 
pears to be accelerating [l], Q has provoked discussion of 
a number of models for extended gravity involving non- 
linear interactions in the Ricci scalar R: 



S 



1 

2k 



d A Xy/=g[R + f(R)} + S n 



(1) 



Here k = 8ttG, S m is the usual matter action and f(R) 
describes the new physics in the gravity sector; setting 
f(R) = — 2 A corresponds to the canonical Einstein- 
Hilbert action in General Relativity (GR) with a cos- 
mological constant A. The idea is that if cosmologi- 
cal data could be fitted by the use of some nontrivial 
function f(R), one might avoid the theoretical difficul- 
ties and fine-tuning issues related to a pure cosmolog- 
ical constant. However, it has been shown that when 
understood as a metric theory, the action |T]) can lead 
to predictions that are not consistent with Solar Sys- 
tem measurements 0, 3 HI ■ While observations require 
a parameter |7ppn — 1| ^ 10 -4 [|| in the Parametrized 
Post-Newtonian (PPN) formalism, the value predicted in 
metric f(R) theories is typically 7ppn ~ 1/2. This is cer- 
tainly the case 0, 0, for the first simple f(R) models 
suggested in the literature @, [1| ■ It is however difficult to 
make a completely generic prediction of this result and 
there have been many arguments both for 0, 0, [f| Q and 
against [H EH, El E3, EL El metric f(R) gravity fail- 
ing Solar System tests. In particular, more complicated 
f(R) functions have since been suggested which claim to 
yield 7P PN « 1 [11 El El El- 

In this paper we set up the conditions which the func- 
tion f(R) must fulfill, so that a solution to the field equa- 
tions which is compatible with Solar System observations 
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exists, in particular with 7ppn ~ 1- However, we will also 
argue that the mere existence of such a solution does not 
imply that a model is consistent with observations. Since 
metric f(R) gravity is a fourth-order theory, spacetime 
geometry and matter are not in as strict a correspondence 
as in General Relativity; depending on the boundary con- 
ditions on the metric, a given matter distribution can be 
consistent with different static spacetimes and with dif- 
ferent values of 7ppn- Moreover, no physical principle 
tells us that only the boundary conditions corresponding 
to 7ppn ~ 1 solutions should be acceptable. The ques- 
tion is then, which solutions are the most natural ones? 
How plausible is it that the collapse of a protostcllar dust 
cloud leads to the formation of the spacetime observed in 
the Solar System? To answer these questions one would 
ideally like to study the full dynamical collapse, and given 
a domain of reasonable initial conditions, determine the 
attractor in the configuration space of possible solutions. 
This computation is beyond the scope of this paper how- 
ever. We will instead approach the problem by study- 
ing how the time stability argument constrains the phase 
space of configurations with the desired properties. 

The conditions that a generic metric f(R) model 
should satisfy in order to yield acceptable solutions are: 
first, the Ricci scalar should closely follow the trace of 
the energy-momentum tensor inside a changing matter 
distribution, where at the same time the dimensionless 
quantities f/R and F = df/dR should remain much 
smaller than 1 at regions of high density. Second, the 
effective mass term m 2 R for a perturbation in the Ricci 
scalar should be positive in order to assure that the GR- 
like 7ppn ~ 1 configurations are stable in time. Third, 
the mass m 2 R should remain small so that a finite domain 
of static, GR-like configurations exist. This is guaran- 
teed if m 2 R < where r Q is the radius of the Sun. 
If this last condition is not fulfilled, the domain of GR- 
like configurations shrinks to essentially a point in the 
phase space, while a continuum of equally good, but ob- 
servationally excluded, solutions still exists. In such a 
case the credibility of the theory requires an argument 
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as to why the particular GR-like configuration should be 
preferred. None of the models so far proposed in the lit- 
erature, including Refs. [H, El, EH , satisfy all of these 
constraints, and we also failed to construct a model that 
would. Largely this failure comes from the difficulty to 
keep both the function F and m 2 R ~ 1/(3^,^) small si- 
multaneously when the Ricci scalar follows the matter 
distribution, R ps tip. 

It should be noted that the above considerations only 
apply for the desired GR-like 7ppn ~ 1 configurations 
when a model is tuned to mimic a (very small) cosmo- 
logical constant. Metric models with a true cosmological 
constant plus some additional sufficiently small f(R) cor- 
rection can be perfectly fine. Thus the above arguments 
do not exclude generic f(R) modifications, such as might 
arise from quantum corrections, to the Einstcin-Hilbcrt 
action. One should also note that it is in general easy to 
construct stable attractor solutions yielding 7ppn ~ 1/2 
in metric f(R) theory. It is the precision data from the 
Solar System which makes these solutions unacceptable. 

The paper is organized as follows. We start by review- 
ing the Solar System constraints in S ec, [Hi We consider 
the Dolgov-Kawasaki time instability [17| in Sec. IHII and 
discuss the corresponding stability criterion for spher- 
ically symmetric configurations. Section IIVI considers 
static configurations and the possibility for metric f(R) 
gravity to follow stable, GR-like solutions that are com- 
patible with Solar System constraints. We find that the 
condition for finding a finite domain of boundary condi- 
tions giving rise to a GR-like metric is nearly orthogonal 
to the time stability condition. Finally, Sec. [V] contains 
our conclusions and discussion. 



II. SOLAR SYSTEM CONSTRAINTS AND THE 
PALATINI TRACK 

Let us begin by reviewing the main constraints from 
the Solar System observations on static solutions in met- 
ric f(R) gravity. Varying the action (Q]) with respect to 
the metric gives the equation of motion: 

{l + F)R pv -\{R + f) gtlv 



(2) 



where F = f jR = df/dR and □ = ff^V^V,,. Taking the 
trace of this equation one finds: 



DF-^R(1-F) 



(3) 



If F — > and / — > R this equation reduces to the stan- 
dard algebraic GR relation between the Ricci scalar and 
the trace of the energy- momentum tensor T . In a generic 
metric f(R) theory R is a dynamical variable however, 
and the theory may exhibit an instability which we will 
discuss in the next section. Assuming a static, spherically 
symmetric metric g^ u , 

ds 2 = g^xW = -e A ^dt 2 + e B ^dr 2 + r 2 dQ 2 , (4) 



the full field equations (|2|) reduce to the following source 
equations for the metric functions A and B in the weak 
field limit (to first order in small quantities): 
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(5) 
. (6) 



where 7 = rF' /2(1 + F), a prime refers to a derivative 
with respect to r, and we have neglected pressure so that 
T ~ —p. The parameter 7 and the terms in the square 
brackets highlight the deviation from General Relativity. 
The value of 7ppn ~ —B/A far away from a gravita- 
tional source depends on the continuous evolution of A 
and B throughout the Sun. It is particularly sensitive 
to the evolution through the core where the density is 
the highest. Hence, to obtain 7ppn ~ 1 and the correct 
gravitational strength in the Solar System, the only so- 
lution, not obviously dependent on an enormous amount 
of fine-tuning [22j|, is that the extra terms in Eqns. (J5JH]) 
must remain small throughout the interior of the Sun. 
Now, if the extra terms can be neglected in the B' equa- 
tion, one finds that B < 10 -6 throughout the interior of 
the Sun [f|. It then becomes clear from the A' equation 
that f/R, F, and rF' need to be very small compared 
to 1. However, to make the correction vanish in the B' 
equation one in addition needs to require that the Ricci 
scalar traces the matter density R/np w 1. So, barring 
perhaps some fantastic fine-tunings, the only possibility 
is that one finds a configuration for which: 



F <C 1 , //i?<l, and Rtzup. 



(7) 



Note that this limit was also discussed in Ref. [12| . Here 
we see that the above conditions are a necessary require- 
ment for fulfilling the local gravity constraints. 

In GR the equation R = —kT is of course exact (re- 
member that we are neglecting pressure throughout so 
that T f=a — p), but this is in general very difficult to ar- 
range in a metric f(R) model. The problem lies in the 
dynamical nature of the Ricci scalar in metric f(R) grav- 
ity. To sec this, consider the static trace equation ([3]) in 
the weak field limit: 



F" + -F' 



(R-Kp- FR + 2/) 
- K p) , 



(8) 



where we have again assumed that pressure is negligible. 
Assume now that F = Fo at the center of the Sun. If the 
nonlinear term ~S(F) is small compared to Kp, then the 
solution for F becomes: 



F(r) 



dr 



,2Gm{r') 
3r' 2 



+ F 



(9) 



where m{r) = dr' '4irr' 2 p. In the case of the Sun this 
implies that F evolves only little: \F(r) - F \ < 10" 6 0]. 
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The solution (|9J) is in general not compatible with Rfanp 
as required by the conditions |(7J). 

Let us now set S(-Fo) = K Po at r = 0. If the gradients 
somehow remain small throughout the evolution, then 
the solution follows the Palatini trace equation: 



£(F) = R(l - F) + 2/ = 8irGp . 



(10) 



For small F and f / R this evolution would be consistent 
with the condition R m up and, since we are following 
the Palatini track, give 7ppn ~ 1 (see e.g. [Hj]). Whether 
such a solution actually exists is more difficult to prove. 
However, one can study under which conditions such a 
solution, if it exists, would be an attractor and whether 
it would also be sufficiently stable in time. 

Finally, let us note that if E(F) <C «p so that F is 
given by Eqn. @, small F and f/R result in a different 
class of solutions with R/ up <C 1. In this case the field 
Eqns. dS][Hl) reduce to 



(rB)> 



D 



A' k 2F' . 

r 



(11) 
(12) 



It is easy to show that together with Eqn. (J9j) these give 
7ppn ~ 1/2. 



III. TIME STABILITY 

Let us first consider the time stability of spherically 
symmetric configurations in generic f(R) models. Per- 
turbing around some arbitrary configuration, R(r) — > 
R(r,t) = R(r) + SR(r,t), and expanding to first order 
in <5_R, SR' and SR, where the prime refers to a derivative 
with respect to r and the dot to a derivative with respect 
to t, one can write the trace equation © in the following 
form in the weak field limit: 



(d?-V 2 )5R = -m%5R + 2 



R 

V 2 i? 



3F 



Fr 



-A 



R'SR' 



F 



,R 



,RR 
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-kT - R(l —F) — 2f and 
1 



'R 



3F 



.(1-F-e), 



{R')\ (13) 



(14) 




FIG. 1: The gradient term proportional to SR' in Eqn. (I16|l . 
normalized to m 2 R /rQ, for various f(R) models (R — up): 
—p 4 /R (solid blue), —p 4 /R + aR 2 /p 2 (dashed green), Hu 
& Sawicki (dot-dashed red) [13 ]. and aR\og(R/p 2 ) (solid 
black). The actual density profile used in all figures corre- 
sponds to the known density profile of the Sun with a central 
density of 150 g/cm 3 and with a roughly exponential depen- 
dence on r. We have also superimposed a constant dark mat- 
ter distribution on the profile of the Sun with pom = 0.3 
GeV/cm 3 . 




FIG. 2: The parameter F given as a function of the radius for 
various f(R) models (R = up): —p 4 /R (solid blue), ~p 4 /R + 
aR 2 1 p 2 (dashed green), Hu & Sawicki (dot-dashed red) [T3 ] . 
and aRlog(R/p 2 ) (solid black). 



background value R(r). Assuming that the configura- 
tion R(r) we are perturbing around is a solution to the 
static equation, the second line in Eqn. Q13p drops out. 
Moreover, for most cases the gradient term proportional 
to 5R' is completely negligible inside a stellar object and 
can be dropped as well. See Fig. [T] for some examples. 
The behavior of the perturbation around a static, spheri- 
cally symmetric solution is thus governed by the equation 



with 



e SR = R(F — F) - 2(f - f) 
'Frr 



+ 3F> fl 



F 



,RR 



F\r 

Fr 



'\2 



(iO 



(15) 



Here a tilde is used to denote that a quantity is per- 
turbed, i.e. it is a function of R(r,t) as opposed to the 



,5R. 



(16) 



Note that the mass m 2 R only depends on the background 
value of the Ricci scalar R(r). 

Table [I] and [Til show the components of m 2 R in some 
particular models and the corresponding parameter val- 
ues used in all figures are displayed in Table IIIIl As 
discussed in the previous section, F needs to be small 
compared to one for GR-like configurations (see Fig. [21 
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we will discuss this constraint further below) . When this 
is the case one typically finds that also e <C 1, so that 
m% « l/3F fl . It then follows that if F >R < 0, then 
mr R < and the coefficient of SR is negative for the 
configuration in question, so that system exhibits an in- 
stability. This is the instability first found by Dolgov and 
Kawasaki in the context of an f(R) = —/i 4 /R model [Tt} 
(see Ref. [IH for a more general case). It is important to 
note that the instability depends not only on the model, 
but also on the particular configuration. Certain configu- 
rations in a given model are more stable than others and 
the instability may even vanish in some cases. 

The nature of the instability is most transparent in the 
special case with constant curvature. Then m R is a con- 
stant and one can obtain an exact solution for 8R(r,t). 
Expanding SR in Fourier modes, one finds that a mode 
with wave vector k has the time dependence 

5R k (k, t) ~ e ±*V k2 + m2 n\ (17) 
so that for negative m R ~ 1/3F, R , all modes with k < 



\m R \ are unstable. This does not necessarily rule out 
a model however. If for example \m R \ ~ H , then the 
instability time is much longer than the lifetime of the 
Solar System and the model is safe. Moreover, whenever 
|mfl | _1 is much larger than the size of the physical system 
under consideration, only modes corresponding to scales 
much larger than the system are unstable and this can 
not alter its local geometry. 

Now, assume that we have a GR-like solution, such 
that R ~ — kT w 8irGp. One then has 

10 (^i 2 ) (g/cm 3 ) ' ^ 

where we have used A w 0.73K/? cr it- Hence, for a pure 
— /i 4 /i? model the mass squared is on the order of 

m -~-< 10 " 26s, " 2 (?) 2 (i7^) ! '■ (19) 

This system is violently unstable at normal densities for 



/Cg) 

4 

R 
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"15" + a — 

R n 2 

! ClCR//i 2 )" 

c 2 {R/y?) n + l 
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aR log ■ 
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R? 



T2+ 2aR 



cm 



4R n + 1 (l + l/c 2 R n ) 2 



a(l + log 



R 



l/3F fl 



6(q - 1/R 3 ) 
c 2 2 R n+2 (1 + l/c 2 R n ) 3 



3cm(n + 1) 1 - ^/ C2 R n 
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3a 



TABLE I: The function F and the dominant term of the mass squared m 2 R = (1 — F — e) /3F.r for different f(R) models, where 
R = R/^ 2 . 
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fi 4 R 2 
R n 2 

1 c 2 {R/fi 2 )" + 1 

aR log ^~ 
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2a 2 fl + 2q/i? 2 -3(l + Kf/i?)/fl 3 + 5/fl 5 4a-l/_fi 3 /1T\2 

a-1/5 3 i? 3 (a- 1/ii 3 )^ 



ra + 2 Z' /tT 
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2ci 



" + 2 for, iARV 
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TABLE II: The parameter e for different /(-R) models, where D = (1— ^j/c 2 R)(l + l/c 2 R n ) and a bar indicates that a quantity 

is dimensionless and measured in units of fi, e.g. R = R/jj 2 . In the third model above, originally suggested in Ref. [l2|], we 
have only kept the leading terms in e. 



5 



all scales larger than ~ 10~ 18 m, if p is fixed to account 
for the present accelerating expansion of the Universe. 

As was pointed out by Dick [lj| and later discussed by 
Nojiri & Odintsov [2(J, adding a conformal term aR 2 / p 2 
can stabilize this system; for f(R) = —p 4 /R + aR 2 / p 2 
the previous approximation for the mass reads: 



R3_ 

6p 4 



1 



1 - aR 3 /p ( 



/£ 

Q 



(20) 



where the last step assumes that aR 3 /p 6 > 1. If a ~ 1 



so that one always finds 
■p 2 /a~ (10 18 s)" 2 . 
However, the above stabilization mechanism runs into 
problems with the conditions in Eqn. ((7|). Indeed, for 



this may be true even for R ~ p 
a very small positive mass m 2 R r> 



1 and R ~ up ^ p one has 



F 



R 2 



R 



»1, 



(21) 



so that the configuration would clearly not be GR- 
like [23|. The problem is that changing F modifies the 
effective strength of the gravitational constant G c s = 
G/(l+F), which controls the buildup of the gravitational 
potential inside the star. In fact, for a ~ 1 the effect is 
so strong that it would weaken the gravitational force so 
much as to prohibit the growth of any density contrasts 
much above the critical density. This argument can be 
turned around to a constraint: in order for the function 
F to remain small inside the densest objects we have rea- 
sonably accurate information on, the neutron stars, one 
has to have aKp nnc \/ p 2 <§; 1, where /9 n uci is the nuclear 
density [24| . Since p 2 ~ Kp C rit we find that 



< 



Pcrit 



10" 



(22) 



This is quite a stringent constraint, but it does not rule 
out the model based on the required time stability. In- 
deed, for example with a = 10 -47 one has a(np) 3 /p 6 = 1 



when 



10 pent- For any density higher than this 



value, the system is stable in time with a very large pos- 
itive mass squared given by the formula in Eqn. (|20[) : 
m 2 R ~ p 2 /a ~ (10 -6 s)~ 2 . There is a caveat to this ar- 
gument however, since for these parameters the gradient 
term proportional to 5R' becomes very large inside the 
Sun (see Fig. ((!])) and the simplified equation (|16p can 
no longer be trusted. 

The complete mass squared function for the model 
with a fine-tuned conformal aR 2 term, using the exact 
expression (|14|). is shown in Fig. [3] (dashed green curve). 
The lower panel displays the absolute value \m 2 R \ and the 
upper panel the sign of m 2 R . In the model at hand the 
mass would remain large and positive throughout the en- 
tire interior of the Sun, which is the necessary condition 
for time stability. The GR-likc configuration does be- 
come unstable at low densities, but this would not nec- 
essarily change the value of 7ppn in the Solar System. 

Fig. [3] also displays m 2 R for several other models (for 
parameter values used in each model see Table ITO|) : the 




FIG. 3: The mass squared m 2 R in units l/r% for various 
f(R) models (R = up): -p 4 /R (solid blue), ~p 4 /R + 
aR 2 /p 2 (dashed green), Hu & Sawicki (dot-dashed red) [T3 |. 
and aRlog(R/p 2 ) (solid black). The horizontal solid gray 
line corresponds to the limit l/r§>. For the —p 4 /R and 
aRlog(R/p 2 ) models we have also plotted 1/3F_r with dot- 
ted blue and dotted black lines, respectively. The upper panel 
displays the corresponding sign of the mass squared where we 
have excluded the —p 4 /R model for which mfj is strictly neg- 
ative. 
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Parameter values 



= 4A/V3 



4A/V3, a = l(T 



p 2 = (8315 Mpc) -1 , n = l 
ci/c 2 = 6 x 0.76/0.24 , 

ci/(£ = 10" 6 x 41 n+1 /n 



4Ae 



(l-a)/a 



a = 1/loglO 3 



TABLE III: Chosen parameter values for the different f{R) 
models in Figs. IH3I where A = 0.73Kp cr it- For the third 
model, originally suggested in Ref. we have used values 
given in the original publication. Note that a value n = 4, 
which was also discussed in [T2| , would result in an even larger 
value of m 2 R in this scenario. 



solid blue line represents the simple —p 4 /R model, which 
has a very large negative mass inside the Sun, and the 
dash-dotted red curve shows the mass function in a model 
by Hu & Sawicki (HS) [3. The HS model fulfills the 
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conditions J7J) by construction, and its very large posi- 
tive mass guarantees time stability. The fact that m R 
becomes negative around r ~ 6r Q in the HS model is 
caused by the e term in the complete expression (|T4"|). 
but this does not necessarily have any effect on 7ppn- 
Moreover, this feature is sensitive to the particular form 
of the exterior density profile (where we have neglected 
for example the Solar wind) and it is not important for 
our main results. Overall one sees that the expression 
m 2 R ~ 1/3Fr is a very good approximation for the first 
three models described in Table U except for a small re- 
gion around r ~ 6r Q where the e term may come into 
play. For the logarithmic model this approximation is 
only good for very low densities and we will discuss this 
in more detail in section Hvl 

We can summarize this section as follows: the stability 
of a static, spherically symmetric GR-like configuration 
with R ks up is predominantly governed by the mass term 
m 2 R ~ I/ZFr. If F R < 0, all perturbations with wave- 
lengths larger than l/|m#| will be unstable so that for a 
large mass the curvature inside a stellar object will 

evolve rapidly before some nonperturbative effect stabi- 
lizes the system. Hence, in order for a model to be stable, 
m R should be positive throughout the Sun for GR-likc 
configurations |25|]. Both the model with a fine-tuned 
conformal aR 2 term (apart from the caveat mentioned 
above) and the HS scenario do satisfy all constraints dis- 
cussed so far. 

IV. AN UPPER BOUND ON m| ? 

As mentioned in the introduction, a given matter dis- 
tribution can be consistent with many different static ge- 
ometries, depending on how the boundary conditions are 
defined at the center of the star. One always requires 
that the exterior metric is asymptotically flat and so dif- 
ferent solutions are characterized by different values of 
7ppn- There appears to be no a priori preference of one 
solution to another and indeed the question is: what is 
the most probable configuration to arise through gravita- 
tional collapse? Lacking a dynamical calculation we are 
restricted here to study how special the GR-like solutions 
are in the phase space. 

Consider a static solution R(r) = Rt(t) + S(r) where 
Rt is the solution to the Palatini trace equation and 
5/Rt <C 1 so that R remains very close to the Palatini 
track. Note that the function S(r) is not a true perturba- 
tion since Rt (r) is not a solution to the complete metric 
trace equation ([3]). However, one can easily obtain the 
equation governing S via Eqn. (jf 3[) . giving 

V 2 <5 = _L(l-F- £ )*-2^i^' 
o-F.R r.R 

-V 2 Rt-^(R' t ) 2 , (23) 

where F and its derivatives are functions of the "back- 
ground" value Rx(r). Similarly, the "perturbed" quan- 



tities in the definition for e are functions of R(r) = 
R T (r)+6{r). 

In analogy with the above analysis for SR(r,t), let us 
consider a constant density object so that Rt = const., 
giving 

V 2 8 = m%6. (24) 

The mass term m 2 R is of course the same mass that ap- 
pears in the equation for SR(r,t). Now, for m 2 R < 0, the 
solution for 5(r) is decaying so that the Palatini track 
acts as an attractor for the solution R(r). This is exactly 
the behavior that was demonstrated by solving the full 
field equations in Ref. for a f{R) = —p 4 /R model. 
However, if m 2 R > 0, the solution for S(r) will also con- 
tain a growing component: 

r r 

The fine-tuning problem we have to face is manifest from 
this equation: the time stability argument of the preced- 
ing section requires that m 2 R > 0, so that Eqn. (|2"5|) with 
its growing mode is what describes the deviation away 
from GR-like solutions. 

From Sec. Mil we already know that setting up a GR- 
like configuration requires fine-tuning at the center of the 
star. However, what Eqn. (|25p implies is worse: starting 
from r = we could always choose a boundary condition 
(Fo,Fq) that kills the growing mode, but as any per- 
turbation around this solution would be exponentially 
enhanced, the boundary condition must be set with an 
incredible precision when m R is large. Numerically such 
a solution can always be found by use of a differentia- 
tion method that kills the growing mode as was done in 
Ref. [l2| . However, different choices of boundary con- 
ditions would lead to other physically, but not obser- 
vationally, acceptable spacctimcs. For example, if one 
starts from a point a little off from the GR track, S ini- 
tially grows exponentially pulling the solution away from 
R « up. Then the nonlinear terms typically become neg- 
ligible in Eqn. and R(F) starts to approach the solu- 
tion of Eqn. For R/np <C 1 this limit corresponds to 
the evolution of A and B given by Eqns. (jlltil 2p . which 
leads to 7ppn ~ 1/2. Thus, for a large m 2 R the nearly 
singular static GR-like solution is surrounded by a contin- 
uum of equally acceptable configurations, however with 
observationally excluded values for 7ppn • 

Hence, given no physical reason to prefer a given set 
of boundary conditions, it would appear more natural 
to expect that the metric around a generic star would 
correspond to 7ppn different from one. As stated above, 
to make a definitive statement would require solving the 
dynamical problem of collapse, but this is beyond the 
scope of the present work. Nevertheless we believe that 
we have identified a potential problem for metric f{R) 
gravity models: for an f(R) model to be credible, it is 
not sufficient to provide a mere proof of existence of a 
GR-likc solution, but one should also give an argument 
as to why this particular solution is preferred. 
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The situation would be ameliorated if the growing 
mode is not excluded, but the length scale dictating 
the growth of the perturbations, 1/niR, is small enough. 
Roughly one should have 

m\ < ^ , (26) 

throughout the Sun. However, both the HS scenario and 
the fine-funded f(R) = -p A / R + aR 2 / p 2 model fail this 
constraint by a large margin, as can be seen from Fig. [3] 
This is also the case for the model in Ref. [ll[ where a sta- 
bilizing conformal term creates a behavior very similar to 
the aR 2 model discussed here. The same argument also 
applies to more recent models introduced in Refs. 
These scenarios behave very similar to the HS model at 
late times, but were designed to also account for infla- 
tion at very high energies. For example, for the model 
suggested in Ref. [T3 ]. 

a R m + l - SR n 
f{R) = l + 7 # ' (2?) 

where the authors chose m = I = n for simplicity, and 
n > 2, one can show that the mass squared is given 
byH 

™«~+3^T) (§)""'■ < 28 > 

Here Rj ~ (10 15 GeV) 2 is set to the scale of inflation, 
and so this mass is enormous in comparison with the 
bound (|26|) inside the Sun. 

The generic problem is that a small value of m 2 R ~ 
1/3-Fij requires a large value for F R , However, at the 
same time one also needs F <C 1 in order to obtain a 
reasonable gravitational potential. This tension is what 
makes it difficult to find a suitable function f(R). Let 
us illustrate the problem further by trying to construct 
an explicit model by rst making sure that the toughest 
requirement is satisfied. At the center of Sun R ss up ~ 
10 31 A ~ 10 -3 /rQ, which is much smaller than the upper 
limit on m 2 R . Thus we can take 

m 2 R ~+R. (29) 

Using this together with the formula F^ R ~ l/m R , we 
can construct a candidate model: 

f(R)=aRlog^, (30) 

where p 2 = 4Ae( 1 ~")/ Q in order to obtain the desired 
accelerating expansion of the Universe at present. Fur- 
thermore, demanding that F <C 1 in the interior of the 
Sun yields a < 0.01 so that p 2 > e 100 A H. So, curi- 
ously enough the Solar System constraints would force 
this model to create the desired accelerating expansion 
without an extremely small energy scale ~ \/A 



Unfortunately, there is a flaw in the above argumen- 
tation, since we implicitly assumed that e <C 1 by as- 
suming m 2 R ~ I/ZF^r. This assumption was fine for the 
discussion in the previous sections, but it fails here. In- 
deed, when F tR is large, the gradient term proportional to 
[3F^ R ](R' / R) 2 in e may also become large (see Table |H|) . 
We can estimate the size of this term using an exponen- 
tial density profile for the Sun p ~ po/(l + e^ r-r °- ) ), 
where £ ~ lOrQ 1 : 

Now, since 3F R > r| one gets [3F R ](R' / R) 2 > 1 and 
our simple estimate for the mass fails. A more careful es- 
timate in the model (|30| finds that the mass is dominated 
by the gradient term and one has 

<«^l gr ad = -(^) 2 ~- 100 4- (32) 

This mass squared is actually negative, so the fine-tuning 
problem is no longer an issue. However, we have recre- 
ated a time instability corresponding to the characteristic 
length scale £ -1 of the system. This behavior is clearly 
visible in Fig. [3] where we have plotted both the full mass 
squared m 2 R (solid black) and the bare function 1/3F^ 
(dotted black) for the aR\og(R/ p 2 ) model. 

V. SUMMARY AND DISCUSSION 

We have shown in this paper that attempts to find sta- 
ble static solutions with 7ppn ~ 1 in metric f(R) models, 
designed to also account for the accelerating expansion 
of the Universe, lead to a string of constraints on the 
model parameters. One must find a configuration for 
which simultaneously F = df jdR and / /R remain small 
compared to one in the interior of the star, where the 
strength of the gravitational field is built up, while the 
Ricci scalar traces the matter distribution: R a up. (See 
also Ref. (HI)- In addition, for the configuration to be 
stable in time, the effective mass term m 2 R for a perturba- 
tion in the Ricci scalar needs to be either positive , or 
if negative, must be much larger than the size of 

the physical system under consideration. Furthermore, 
we showed that unless m 2 R < 1/Vq, the domain of GR- 
like static configurations shrinks to essentially a point in 
the phase space, while for example a continuum of solu- 
tions corresponding to 7ppn ~ 1/2 exists. 

Hence, in particular for large positive m R , it would 
appear more natural to expect that the metric around 
a generic star would correspond to 7ppn different from 
one. To make a more definitive statement one should 
solve the dynamical gravitational collapse of a protostel- 
lar dust cloud, which is beyond the scope of this paper 
however. Moreover, to a degree the problem with the 
boundary conditions would merely be translated to set- 
ting the initial conditions for the collapse. Nevertheless, 
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we believe that we have identified a potential problem 
in that to make a given metric f(R) model credible, one 
should give an argument as to why the GR-like configu- 
rations should be preferred. Otherwise, if the PPN and 
stability constraints are supplemented by our fine-tuning 
argument, it seems unlikely that any f(R) model can 
pass the test - unless one gives up the hope that the the- 
ory is also responsible for the accelerating expansion of 
the Universe. This is because the condition rri 2 R < l/Vg, 
combined with m 2 R ~ 1/3-Fjj, implies that Fr needs to 
be large. However, since at the same time one needs 
F « 1, a tension is created that makes finding a suitable 
function f(R) very difficult. 

Let us finally note that while completely GR-like con- 
figurations are hard to construct in f(R) models, it does 
not mean that such theories would be somehow funda- 
mentally ill. In fact many metric f(R) theories could de- 
scribe gravitational physics quite well in most situations; 
it is the very precise information on 7ppn from Solar 
System experiments which eventually forces one to set 



R ?s up. Indeed, if one looks even at the simplest model 
with f(R) = —p 4 /R, one finds that setting R essentially 
to any other value than up gives 7ppn ~ 1/2 More- 
over, whenever one has R ~ p 2 , one finds m R ~ p 2 , so 
that these configurations are effectively free of any stabil- 
ity problems. Thus, one sees that the Dolgov-Kawasaki 
instability and fine-tuning problems depend not only on 
the theory, but also on the particular given configura- 
tion within a given model. Even in the simplest —p^/R 
model one can construct sufficiently stable spacetimes for 
stellar objects that are consistent with the accelerating 
expansion of the Universe; these are only excluded from 
describing reality by the very precise PPN constraints. 
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By this we mean that it is in principle possible that an 
evolution of B/A quite different from the one in GR, 
could actually lead to the same value for the exterior 
of a star. However, away from the GR track the outcome 
becomes sensitive to the form of the density profile, which 



9 



leads to even more uncertainties as to how the actual 
dynamical gravitational collapse would proceed. 

[23] In general, F > 1 implies that the leading approximation 
in Eqn. (|20p is no longer valid. However, it still holds for 
f(R) = —fx 4 /R + aR 2 /fx 2 , since the contribution of the 
conformal term in F will cancel with the leading term in 
e where additional terms remain small compared to one. 

[24] The loophole to this argument is if one only follows R ~ 
up inside the Sun, but that the Ricci scalar is allowed 
to significantly deviate from this relation inside neutron 
stars. In such a case one could maintain F -C 1 for much 
larger values of a, resulting in a smaller m 2 R for the Sun. 



This scenario seems very contrived however. 
[25] Configurations with a small enough negative m 2 R can be 
accepted as well. The minimum requirement is then that 
at least all perturbations with wavelengths smaller than 
the size of the physical system under consideration should 
remain stable. 

[26] However, note that since F — a(log(R/ fj, 2 ) + 1) = 
a log (i?/4A) — 1 + 2a, there is some tension between get- 
ting F < 1 inside both the Sun and a neutron star. Al- 
though this may indeed prove relevant in further con- 
siderations, it is not of importance for the discussion at 
hand. 



